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The calculations of the light curves of thermonuclear supernovae are carried out by a method 
of multi-group radiation hydrodynamics. The effects of spectral lines and expansion opacity are 



taken into account. The predictions for UBVI fluxes are given. The values of rise time for B and 
V bands found in our calculations are in good agreement with the observed values. We explain 
why our results for the rise time have more solid physical justification than those obtained 
by other authors. It is shown that small variations in the chemical composition of the ejecta, 
produced in the explosions with different regimes of nuclear burning, can influence drastically 
the light curve decline in the B band and, to a lesser extent, in the Fband. We argue that recent 
results on positive cosmological constant A, found from the high redshift supernova observations, 
could be wrong in the case of possible variations of the preferred mode of nuclear burning in the 
earlier Universe. 



Introduction 



For decades, Type la supernovae (SNe la) are treated as very suitable objects for distance 
measurements in the Universe and for deriving cosmological parameters (Sandage, Tammann, 
1997; Ruiz-Lapuente, 1997). There are several reasons for this. The first one is the intrinsic 
brightness of the SNe la. Due to their luminosity, one can measure an appreciable flux even 
if they are at high redshift, like z ~ 0.5 — 1. Besides that, the spectra and light curves of 
SNe la seem very uniform at first glance. However, more thorough investigation demonstrates 
the diversity within this class of objects (Pskovskii 1977, Bartunov & Tsvetkov 1986, Branch 
1987, Barbon et al. 1989, Phillips et al. 1987; Filippenko, et al. 1992, Leibundgut, et al. 1993). 

Quite a while ago, Pskovskii (1977) has shown that there is a correlation between the maxi- 
mum luminosities of nearby SNe I and their post-maximum decline rates. This dependence was 
confirmed subsequently from the observations of many low-z supernovae and studied by several 
workers (Bartunov & Tsvetkov 1986, Phillips 1993, Hamuy et al. 1995,1996a, Hamuy & Pinto 
1999). 

Modern observational technique allows discovering and studying distant, high-z, supernovae. 
The first results were obtained by N0rgaard-Nielsen et al. (1989): a few years of their observa- 
tions have brought up only two supernovae at z ~ 0.3. 

Presently, there are several search groups in the world which observe distant supernovae with 
the largest terrestrial telescopes (and sometimes with the help of the Hubble Space Telescope). 
The observational technique is advanced to the level allowing the discovery of dozens high- 
z supernovae per one observational period of a couple of weeks. The observational material 
produced by those groups allowed to estimate the cosmological parameters such as the Hubble 
constant Hq, the matter density £l M and the vacuum energy Q\ in units of the critical density 
(for definitions see, e.g., Carroll et al. 1992). The dependent quantities, like the deceleration 
parameter q, the ratio of the local value of Hq to the global one, etc., can be derived based 
on this. E.g., Kim et al. (1997) have estimated Hq using the first seven SNe la at z > 0.35 
and refuted the suspicion that the local value of Hq is appreciably larger than the average one. 
Perlmutter et al. (1997) have estimated a probable ratio of the Om and Oa- 

Quite recently, several interesting papers appeared based on a richer statistics of distant 
supernovae (Schmidt et al. 1998; Garnavich et al. 1998a, b; Riess et al. 1998; Perlmutter et al. 
1999). An extremely important result is claimed in that work: the analysis of the observations 
implies with high confidence that the expansion of the Universe is accelerating at present epoch. 

One should note, though, that those results on distant supernovae were based on the maxi- 
mum luminosity - decline rate relation derived from the analysis of nearby objects (in the work 
by Perlmutter et al. 1999, the high redshift SNe la are also used in obtaining such a relation, 
but this does not help in excluding possible effects of evolution, see Drell et al. 1999). Even for 
the nearby SNe la the deviations of individual objects from the relation cannot be explained 
solely by observational errors. The physical understanding of the peak luminosity - decline rate 
correlation is crucial for estimating the cosmological results obtained with SNe la. 

In principle, the slower decline rate for the higher maximum luminosity can be explained 
from theoretical point of view: both are caused mainly by 56 Ni produced during the explosion. 
Radioactive decay of 56 Ni forms the SN la light curve and determines its maximum luminosity. 
On the other hand, large amount of nickel should enhance the opacity of matter. It takes 
longer time for radiation to diffuse through the stellar matter and the light curve becomes less 
steep. However, the decline rate is influenced not only by the amount of nickel (and of other 



heavy elements), but also by their distribution inside the expanding star and by the velocity 
distribution of ejected gas. These distributions are dependent, in turn, on the mode of burning 
propagation in the star. 

After the pioneering work by Arnett (1969), Ivanova et al. (1974), Nomoto et al. (1976) 
the theory of thermonuclear burning in supernovae has been advanced appreciably, yet many 
problems are still not solved, see, e.g., Niemeyer and Woosley (1997), Reinecke et al. (1999), 
Niemeyer (1999). A plenty of models for SN la explosions were built, with various masses (Chan- 
drasekhar and sub-Chandrasekhar) , burning modes (detonation, deflagration, and a variety of 
their combinations) in a range of explosion energies and expansion velocities. As a result of 
dissimilar burning, chemical elements are produced in various proportions and distributed dif- 
ferently throughout the ejecta, even if conditions in a star prior to the explosion are identical. 
This leads to noticeable variations in theoretical light curves. From the comparison of calculated 
light curves with the observed ones one can judge which mode of explosion is preferred in na- 
ture. Different explosion scenarios imply a scatter in the dependencies between the observable 
parameters of the burst, and could explain the events deviating from the standard relations to 
a larger extent than a margin anticipated from the observational errors. 

An extensive set of the models for light curves of SNe la was produced by Hoflich et al. 
(1996,1997), giving some theoretical insight into the observed correlations of the peak luminos- 
ity and the post-maximum decline rate. We have undertaken an independent study of Chan- 
drasekhar and sub-Chandrasekhar SN la light curve models. The results will be published 
elsewhere (Sorokina et al. 1999), some of our calculations for sub-Chandrasekhar models are 
discussed by Ruiz-Lapuente et al. (1999). In this report we present the results of the calcula- 
tions of the light curves for two well-known SN la models: the deflagration model W7 (Nomoto 
et al. 1984) and the delayed detonation one, DD4 (Woosley & Weaver 1994). They both are 
Chandrasekhar mass models and contain similar amount of 56 Ni, but the latter is distributed 
differently due to different burning modes. Our result is that the two light curves are not similar 
in their decline rate, despite the close absolute luminosities at maximum light. The S-band flux 
of W7 goes down too slowly and agrees with the observations of typical SNe la only marginally. 
We explain why the results of other authors somewhat disagree with ours for this model. 

Discussing the possibility of the use of SNe la in cosmology, we suggest that the available 
statistics of distant supernovae does not yet allow drawing firm conclusions on the geometry of 
the Universe. In terrestrial experiments, one cannot predict with certainty the outcome of an 
explosion. The situation for supernovae can be similar: it is quite likely that initial conditions 
influence only the probability of the development of a certain burning mode, but do not determine 
it exactly. Since the light curve shape depends strongly on the regime of burning, one cannot 
predict deterministically the decline rate knowing only the initial conditions. The probability 
of the specific decline rate, which is crucial for obtaining the cosmological parameters, would be 
found only with a sufficiently rich statistics of SNe la at various redshifts. 

Numerical method 



The approach used for light curve modeling in our study is the multi-energy group radiation 
hydrodynamics. The code we have used is called Stella (Blinnikov &: Bartunov 1993; Blinnikov 
et al. 1998). It solves the time-dependent equations for the angular moments of intensity 
averaged over fixed frequency bands, using up to ~ 200 zones for the Lagrangean coordinate 



and up to 100 frequency bins (i.e., energy groups). This number of frequency groups allows one 
to have a reasonably accurate representation of non-equilibrium continuum radiation. There 
is no need to ascribe any temperature to the radiation: the photon energy distribution can be 
quite arbitrary. 

The coupling of multi-group radiative transfer with hydrodynamics means that we can obtain 
UBVRI fluxes in a self-consistent calculation, and that no additional estimates of thermalization 
depth as in the one-energy group models are needed. Variable Eddington factors are computed, 
which fully take into account scattering and redshifts for each frequency group in each mass 
zone. The parameters of the decays of 56 Ni to 56 Co and of 56 Co to 56 Fe are taken from the 
work of Nadyozhin (1994). The positron energy is deposited locally. The gamma-ray transfer 
is calculated using a one-group approximation for the non-local deposition of the energy of 
radioactive nuclei. Here we follow Swartz et al. (1995) and we use purely absorptive opacity in 
the gamma-ray range. Swartz et al. addressed this question using in particular the W7 model. 
They found a very weak dependence on the optical depth and derived the error of only 10% in 
heating rate due to the gamma deposition in the SN la model up to 1200 days after the explosion 
for the effective opacity k 7 = 0.050Y e (Y e is the total number of electrons per baryon). 

Local Thermodynamic Equilibrium (LTE) for ionization and atomic level populations is 
assumed in our modeling. In the equation of state, LTE ionization and recombination are taken 
into account. The effect of line opacity is treated as an expansion opacity according to the 
prescription of Eastman & Pinto (1993). 

The main limitation of our current light curve code is the LTE approximation. To simulate 
some of the non-LTE effects we used the approximation of the absorptive opacity in spectral 
lines, following the results of non-LTE computations by Baron et al. (1996), see the discussion 
in Blinnikov et al. (1998). We do not pretend yet that the light curves in this approximation 
are reproduced absolutely correctly, especially later than ~ 2 months after the explosion. But 
our experience in the light curve modeling of other types of supernovae (Blinnikov et al. 1998, 
Blinnikov 1999, Blinnikov et al. 1999) shows that our results are very reliable for a month or so 
after the maximum light. One can see in Fig. [| below, that the photosphere in B band for the day 
35 in our models is still in the outer layers, so the bulk of the flux is born in the conditions close 
to LTE. What is also important, we are interested here mostly in relative changes of the light 
curves for two classical regimes of burning. We believe that the relative difference is reproduced 
quite reliably in our computations. 

Burning modes in thermonuclear supernovae 

Currently, there is no doubt that the light of SNe la is produced in the decays of 56 Ni to 
56 Co and then to 56 Fe. The sufficient amount of initial 56 Ni as well as the required energy of the 
explosion can be naturally produced during the thermal runaway in a degenerate carbon-oxygen 
mixture. Yet, the regime of the nuclear burning in supernovae is still a controversial issue. 

Arnett (1969) was the first to model the supersonic combustion, i.e. detonation, in super- 
novae. Later, Ivanova et al. (1974) obtained a sub-sonic flame (deflagration) propagating in a 
spontaneous regime with pulsations and switching subsequently to the detonation. Nomoto et al. 
(1976) modeled deflagration propagating due to a convective heat transfer with a parameterized 
flame speed. Both detonation and deflagration have their merits and problems in explaining the 
supernova phenomenon (see, e.g., the review of Woosley & Weaver 1986). 



A prompt detonation, born near the center of a star, would incinerate almost the whole 
star to ashes consisting of iron peak elements (Arnett 1969). This is in clear contradiction with 
observations which show that intermediate mass elements are abundant in the ejecta. Thus at 
some stage (if not always) the burning must be sub-sonic. Here another problem arises: it is 
clear that the combustion must be much faster than suggested by an analysis of propagation of 
a laminar one-dimensional flame (Timmes & Woosley 1992). 

From the microscopic point of view the one-dimensional nuclear flame is a wave described 
essentially in the same way as that formulated by Zeldovich & Frank-Kamenetsky (1938) in 
spite of complications introduced by nuclear kinetics and a very high conductivity of dense 
presupernova matter. It is found that the conductive flame propagates in a presupernova far 
too slowly to explain the supernova outburst correctly: the flame Mach number is of the order 
of a few percent and less (Timmes & Woosley 1992). 

A natural way to accelerate the fuel consumption is the development of instabilities inherent 
to the flame front. As explained in the classical paper by Landau (1944), the hydrodynamic 
instability leads to wrinkling or roughening of the front surface, i.e., to increasing of its area 
with respect to the smooth front, and consequently to acceleration of the flame propagation. In 
extreme cases, the instabilities can lead to a transition from the regime of slow flame propagation 
to the regime of detonation. The role of the Landau instability in supernovae was pointed out 
by Blinnikov et al. (1995) and Niemeyer & Hillebrandt (1995b). As the instability grows, the 
front becomes wrinkled and fractal. The detailed consideration of non-linear stage of Landau 
instability and the calculation of the fractal dimension of the flame front for that case is given 
by Blinnikov & Sasorov (1996). 

Since the flame propagates in the gravity field and the burned ashes have lower density than 
the unburned fuel, the Rayleigh-Taylor (RT) instability is often considered to be the main pro- 
cess governing the corrugation of the front (see Miiller & Arnett 1986, Woosley 1990a, Woosley 
1990b, Livne & Arnett 1993, Khokhlov 1993). The RT instability gives birth to a turbulent 
cascade providing for an acceleration to the flame and additional difficulties in modeling the 
SN event (Niemeyer & Hillebrandt 1995a, Niemeyer 1995, Niemeyer &: Woosley 1997). The 
turbulent mixing of hot products of burning with the cold fuel might trigger a detonation after 
a period of slow expansion of the partly burned degenerate star (Niemeyer & Woosley 1997). 
Here the transition to the detonation may happen via a stage of supersonic spontaneous flame 
propagation (Blinnikov & Khokhlov 1986). However, recent computations by Reinecke et al. 
(1999) show that the RT instability alone is not sufficient for the explosion: it drives a flame 
front too weakly and the needed subsonic speed is not achieved. Additional sources for the 
turbulence or alternative routes for accelerating flames seem to be needed (Niemeyer 1999). 
Even if detonation is established, it cannot be described by a classical Zeldovich-von Neumann 
wave, which is unstable in the conditions of SNe la explosions (see Imshennik et al. 1999 and 
references therein). 

It is clear that a full direct modeling of the burning is not possible because the involved scales 
range from the thickness of a flame ~ 10 -5 cm up to the dimension of a white dwarf, ~ 10 9 cm. 
Moreover, in all cases where turbulence is involved, the chaotic behavior sets in inevitably, and 
a subtle difference in initial conditions can lead to drastically unsimilar outcomes. This is well 
known from terrestrial experiments: after a series of tens or even hundreds cases of slow burning, 
sometime a detonation sets in when all the conditions of the experiment are quite the same (e.g., 
Mader 1979). 



Theoretical Models of Supernovae and Light Curves 



Initial models corresponding to two different regimes of burning were used in our calculations, 
namely, the classical deflagration model W7 (Nomoto et al. 1984) and a delayed detonation 
model DD4 (Woosley & Weaver 1994). The main parameters of these two models are compared 
in Table [j], and their resulting chemical compositions, in Table g. Figures [j], [2] display the 
distributions of density, velocity, and basic chemical elements in the ejecta for those two models. 
It is clear from the figures and tables that the hydrodynamical parameters are approximately 
the same for both models, while differences in chemistry are due to the relevant burning mode 
in a white dwarf (though, the amount of 56 Ni is almost the same in both models, so maxima in 
luminosity are close). Both W7 and DD4 are simple ID models, and each individual SN la can 
be a more complicated 3D fractal event in reality. Yet our computed light curves have a good 
resemblance to the observed ones and reproduce the broad band fluxes quite well. 

The curves in the Fig. || display the numerical results for the models W7 and DD4, while 
asterisks show the observational points for SN 1998bu (Suntzeff et al. 1999; see also Meikle <fc 
Hernandez 1999, and Jha et al. 1999) and crosses the template fluxes for SN 1992A (Hamuy et 
al. 1996b). The results demonstrate clearly that in spite of similar maximum luminosities, the 
models have very unequal decline rates in the B filter, though they are much more similar in 
the Fband. 

It is not easy to explain this scatter in the predicted decline rates by a single factor. The 
outgoing flux of the SN la explosion is influenced by many physical parameters. The differences 
in all parameters for the models are small for each instant. Yet the combination of the minute 
differences seems to be sufficient to change the light curve shapes. For example, the distributions 
of 56 Ni in outer layers, as well as lighter elements closer to the center, are not the same, as shown 
in the Fig. ^. This causes slight changes in gamma-ray deposition and in diffusion of soft photons, 
thus implying unsimilar temperature distributions. The latter are shown in Fig. || for the day 
35 after the explosion. The solid dots mark the layers where the optical depth becomes equal to 
2/3, if measured from the outer boundary of the ejecta. The temperature for W7 is 10% higher, 
and the radius is also a bit larger than for DD4 in this point. This can explain quite well the 
differences in B fluxes for this epoch (see Fig. |3|). 

Let us compare our computed light curves with observations of real supernovae. There 
are several quantitative relations to describe their shapes. One of the most popular is the 
dependence suggested by Phillips (1993), M mSiX (Amf 5 ), where Amf 5 is a drop of the B flux (in 
stellar magnitudes) 15 days after the maximum. The linear regression fit 

M max = a + 6Amf 5 

is usually used for the interpretation of observations. The values of a and b vary strongly, when 
obtained with different sets of SNe la by different authors. Moreover, the dispersion of their 
values is large even within one set. We are more interested in the value of b, the slope of the 
light curve. Its variation is rather large, e.g., Phillips (1993) obtains b B = 2.698 ± 0.359 from 
the observations of 9 nearby supernovae. In contrast, analysis by Hamuy et al. (1995) of only 
six supernovae from the same set yields b = 1.624 ± 0.582. This scatter cannot be explained 
only by the observational errors, they are not so big. See, however, Hamuy et al. (1996a) 
where they obtain much smaller scatter and even lower value of b = 0.784 ± 0.182 using 26 
Calan/Tololo distant supernovae. The value of a is influenced by the Hubble parameter Hq: 
a = a + 51g(flo/85) for Phillips (1993), and a = a + 51g(fT /65) for Hamuy et al. (1996a). In 



reality, supernovae do not lie on a straight line, they occupy an extended region on the diagram 
M max (Amf 5 ). Most probably, this reflects real variations of the explosion physics. One can 
come to the same conclusion after analyzing the correlations between other key parameters of 
SN la, like 56 Ni mass, peak bolometric luminosity, rise time to the maximum light, explosion 
energy, etc (Nadyozhin, 1997). 

The slope of our calculated light curves (especially for DD4 model) and the luminosity are 
comparable with the observable values. We can estimate the decline rate from the results of 
Phillips (1993) and Hamuy et al. (1995) with the formula Amf 5 = (M^ ax — a)/b. Comparing 
this with the rates for W7 and DD4 in our computations we find for Hq = 85 km/s/Mpc: 

B _ j 1.04(Phillips) 0.94(Hamuy95) 1.07(model DD4) for M^ ax = -18.9; 
mi5 ~ | 1.1 2 (Phillips) 1.06(Hamuy95) 0.66(model W7) for M^ ax = -18.7. 

So, we see that the model W7 behaves in the sense opposite to the Pskovskii-Phillips relation: 
being a bit weaker than DD4 it declines slower, instead of fading down faster. 

If we take more recent data from Hamuy et al. (1996a) and Hq = 65 km/s/Mpc, then we get 
higher values for Amf 5 , predicted by 'inverse Phillips' relations, Amf 5 = (M^ ax — a)/b, for the 
M^ ax found in our modeling, but this reflects only the fact that both models are underluminous 
for the long cosmological distance scale. We are interested in relative differences between two 
models, so it is more instructive to compare the absolute magnitudes from the direct formula 
M maj = a + 6Amf 5 . We present the comparison for the two values of the Hubble parameter. 
For H = 85 km/s/Mpc: 

B { -18.8(Phillips) -18.7(Hamuy96) -18.9(model DD4) for Amf 5 = 1.07; 

max ~ | -19.9(Phillips) -19.0(Hamuy96) -18.7(model W7) for Amf 5 = 0.66 

Thus, for the short distance scale DD4 is OK, while W7 will look underluminous. For Hq = 65 
km/s/Mpc: 

B \ -19.4(Phillips) -19.3(Hamuy96) -18.9(model DD4) for Amf 5 = 1.07; 
max ~ j -20.5(Phillips) -19.6(Hamuy96) -18.7(model W7) for Amf 5 = 0.66 

Now, both models are too weak, and this means that for the long distance scale a standard 
model, like DD4, which is 0.4 magnitudes underluminous according to the peak luminosity- 
decline rate relation of Hamuy et al. (1996a), must have perhaps 30 to 40% more 56 Ni (Arnett 
1982, Cappellaro et al. 1997, Mazzali et al. 1998). At the same time W7 is underluminous 
already by 0.9 magnitudes. If we assume the same relative difference between delayed detonation 
and deflagration models with the same enhancement of 56 Ni abundance, then the W7-like events 
will still look 0.5 magnitudes underluminous. This can have drastic implications for the use of 
SNe la in cosmology. 

One can see in the Fig. |3| that the behavior of both calculated light curves is qualitatively 
correct when compared with a couple of well studied individual SNe la. In particular, the rise 
time corresponds to the observed values: the flux of observed SNe la grows up to the maximum 
during 15-20 days (Maza 1981, Pskovskii 1984, Bartunov & Tsvetkov 1997). As it is shown in 
Table || our calculations predict 16-17 days for this value in the B band, and this is in much 
better agreement with the observations than yielded by previous computations of Hoflich et al. 
(1997). (Their fluxes grow for only 10-15 days for the majority of models, including W7.) One 
can also note the realistic behavior of our computed light curves in different spectral bands. The 
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decline rate is less steep for redder bands in both models (see a simple analytical explanation 
for this given by Arnett 1982) But it is also evident that the deflagration model fades down in 
the -B-band too slowly, and does not fit the decline rates for typical type la supernovae, while 
the delayed detonation model agrees with the observations very well. 

Until recently W7 has been in frequent use in many papers where light curves of SNe la 
were presented, as well as when supernova remnants were calculated subsequently. The reason 
for such a popularity is that this model seemed to reproduce the observed fluxes quite well 
(see, e.g., Hoflich 1995 and Eastman 1997). Our results for W7 (especially in the B band) 
differ from the ones obtained in the previously published papers. This can be explained most 
probably by the differences in including the expansion opacity (Karp et al. 1977) into the energy 
equation. The correct treatment of the energy equation is described by Blinnikov (1996, 1997) 
who has shown that there is no effect of expansion opacity when line extinction is taken into 
account in the energy equation for photons. The expansion effect is important only in the flux 
equation. Extinction in lines is much stronger than other sources of opacity for the SNe la. So 
when the effect by Karp et al. is included into the energy equation, this leads to unjustifiable 
enhancement of the energy exchange between radiation and matter. The photon energy is in this 
case artificially redistributed towards longer wavelengths, and escapes the system much easier. 
As a result, the flux from a supernova grows up and falls down faster, and the slope of the light 
curve becomes steeper. 

To check the importance of this effect, we have calculated the model W7 once more, but 
this time we treated the extinction (including the one in lines and the expansion effect) as a 
pure absorption. In such a way, we have made the energy exchange stronger, and so we have 
spoiled the energy equation. Because of this, the decline rate of the light curve in the B band has 
become "better" , although the rise time is shortened and its agreement with the observations is 
much worse. One can notice here that the results of the calculations of our "spoiled" W7 model 
resemble those by Hoflich. Fig. shows how the enhanced energy exchange influences the light 
curve of W7. It is easy to see from the figure that although the flux in shorter wavelengths 
could be called "good" (but with shortened rise time; see also Table ||), it would be absolutely 
untrue regarding the flux in /. The secondary maximum in / band, which is observed in almost 
all real SNe la, has disappeared in the new light curve. Physical processes responsible for the 
presence of this feature in the light curves should yet be investigated. Nevertheless, the lack of 
the secondary maximum clearly demonstrates that some important processes were not included 
(or were included incorrectly) in the treatment of radiative transport. 

Thus, we are able to reproduce the results of other researchers when the wrong energy 
equation is used in our code. When the appropriate physics is included instead, the decline rate 
of the deflagration model W7 is too slow, and does not correspond to the behavior of a typical 
SN la. Two possible reasons for such a result could be proposed: either the real explosions 
are not due to the pure deflagration or, most probably, too many simplifications were supposed 
while constructing W7. More complicated deflagration model could, perhaps, describe the reality 
better (cf. Niemeyer & Woosley 1997, Reinecke et al. 1999). 

Using SNe la in Cosmology 



Due to the maximum light - decline rate dependence, SNe la are now in active use for 
finding cosmological distances. Their observations can help to determine the geometry of the 



Universe, as well as to check the viability of some cosmological models. The idea is to compare 
the photometric distance to a supernova with its value in a given cosmological model when its 
redshift z is measured. The simple formula d p h = {L/AttF) ' defines the photometric distance 
to an observed object if its absolute luminosity L is known. At the same time, the photometric 
distance can be calculated theoretically: if z is fixed, then d^ is a function of the cosmological 
parameters Hq, Qm and f^A (Carroll et al. 1992): 
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Here fi^ = 1 — ^m — ^A, and for $7^ < the hyperbolic sine (sinh) transforms, according to 
elementary formulae, to the trigonometric sine (sin), while \/^fc should read \/|f2jfc|. For Qj. = 
the limit f2fc — ► is easily taken, so sinh disappears from the expression for dth, and only the 
integral is left. 

At present, several search groups in the world do observe the high-z supernovae. The peak 
luminosity - decline rate relation is used by them to estimate the cosmological parameters from 
these observations. The relation is calibrated with the help of nearby supernovae samples, since 
rich data are available in this case (see, e.g., Phillips 1993), but it is used to interpret the 
observations of the much more distant supernovae as well. We already discussed that even 
nearby examples are not fully homogeneous, though they seem to be equivalent at first glance. 
The investigations of the influence of the evolutionary effects on the light curve shape are only 
at the beginning so far. 

The maximum light - decline rate relation could be different for younger galaxies. To un- 
derstand how it could change, one should know how the shape of the light curve depends on the 
initial chemical composition of the exploded white dwarf. And what are the initial conditions 
in general, what are the progenitors of supernovae, do those conditions change with the age of 
Universe? 

The influence of the initial conditions just before the explosion on the emission of the SNe la 
was studied by Dominguez et al. (1999). The role of evolution of SNe la is discussed by Riess 
ct al. (1998), Schmidt et al. (1998) and Drell et al. (1999). (See also the discussion on different 
pre-supernova models in von Hippel et al. 1997, Livio 1999, and Umeda et al. 1999.) 

Schmidt et al. (1998) compare the distances to 8 SNe la in early type galaxies with 19 in late 
type ones and find no significant difference. But here are two points that need be yet clarified. 
First, it is hard to select the progenitors of the same ages in the late type hosts. Second, the 
statistics of such progenitors is too small. One has to build distributions of SNe la at different 
z with good statistics and to compare them. The statistical analysis by Drell et al. (1999) 
shows convincingly that it is virtually impossible to distinguish the effects of evolution from 
the effects of the geometry of the Universe expressed through ^m and 0,\. Drell et al. (1999) 
write that without physical understanding, the additional data, such as source spectra, cannot 
be used to argue against evolution, while observers pretend that there is no compelling evidence 
for evolution in the SNe la samples, based on these data. The physical insight in the problem is 
being developed, e.g. by Hoflich, Wheeler, and Thielemann 1998, Dominguez et al. 1999, Umeda 
et al. 1999. 

Dominguez et al. (1999) find that due to the changes in metallicity and C/0 ratio, the light 
curve shape can also change. This leads to the drift of ~ 0.2 m in distance modulus determined 
with the averaged luminosity - light curve shape relation. This effect was considered for only 
one set of delayed detonation models, and the results were found to be influenced mostly via 



the change of 56 Ni mass. 

All previous work on the evolution of the SN la progenitors is based on the assumption that 
chemistry and hydrodynamic structure of the pre-super novae exactly fix the behavior of the light 
curve. This is a deterministic approach. With this approach the net effect is not necessarily 
large. As we find in this work, the supernova emission may also depend drastically on the mode 
of burning, and this is most probably the factor that cannot be predicted in a deterministic way 
in principle. The main difference in our approach is a conjecture that even for the same initial 
conditions (or minute differences in the initial conditions) the outcome of the explosion can be 
vastly different. For the same chemistry and for the same age of progenitors there can be a wide 
distribution of different SN la events, and this is confirmed by the observations (see, e.g, Hamuy 
et al. 1995, 1996a). If this distribution (there is no sense to speak about individual events!) 
changes with z then one can obtain a wrong value of A. Fig. |3] is a bright illustration of the 
possibility to erroneously determine the absolute flux from a SN la and, therefore, its distance. 
It shows the strong deviation in the decline rate among the SNe la with similar maximum light 
due to the differences in the burning regime, which is hardly predictable in principle, as was 
discussed above. One can only predict a trend for the preferred mode of burning, which can 
indeed depend on the age and the chemistry of the progenitor. 

Let us assume that W7 and DD4 represent two limiting cases for real Chandrasekhar-mass 
explosions. There can be a continuous set of supernovae between those limits, or there can be a 
clustering to the two extreme cases. This does not matter for our purpose now. What matters 
is the indeterminacy principle governing the behavior of the flame inside a star. Thus, when the 
initial conditions are fixed, there can be a trend in the distribution of the real events to be more 
"DD4-like" or more "W7-like". If, for instance, we have more W7-like events at high z then 
we can obtain a wrong result when we try to derive the distance, and hence, the cosmological 
constant from SNe la, based on the local calibration samples, because of different luminosity - 
slope dependence in the past. 

In the recent work on the cosmological supernovae (Riess et al. 1998, Perlmutter et al. 1999) 
the conclusion is reached, that the distance to the high-z supernovae is about 10-15% larger 
than expected for the model with f^M = 0.2, £l\ = 0. It is claimed that a model with f^ > 
is needed to avoid such a deviation. In reality, A could be exactly zero if the fraction of slowly 
burning supernovae would be larger for higher z. The true luminosity would be then smaller than 
found from the peak luminosity - light curve shape calibrations for the local sample, and this 
could mimic the effect of positive A. We are persuaded that so far there is no sufficient statistics 
of distant supernovae to make firm conclusions on the geometry of the Universe. Not only the 
non-zero A can reconcile the theoretical distances dth with the observed photometric ones, but 
also the evolution of the probability of the burning regimes, making the W7-like explosions in 
the present epoch less likely than in the younger Universe. 

Conclusions 



We have chosen as initial configurations for SNe la two models of carbon-oxygen white dwarfs 
with different regimes of burning during the explosion, and simulated the light curves till ~ 2 
months after the burst. We find that the post-maximum decline rate is very sensitive to the 
burning mode. Even minute differences in the 56 Ni distribution seem to be sufficient to heat the 
material differently inside the ejecta and to produce different slopes of the light curve thereby, 



despite the similar maximum brightness. 

The combustion theory on its modern level is not able to predict the mode of burning if the 
initial conditions are fixed. It is conceivable that this prediction is impossible in principle. The 
different burning modes, in turn, can result in unlike temperature distributions. If this is true, 
then the initial conditions do not determine the absolute maximum - decline rate dependence 
exactly. They fix only the probability to obtain one or another slope of the light curve for the 
fixed maximum brightness. This probability can be changed by variation in the initial conditions. 

One can speculate that the future combustion theory will be able to predict only a proportion 
between "fast" and "slow" supernovae, and this proportion can be very sensitive to the initial 
conditions. Astronomers should not wait for such an advancement of the theory, but could 
already now check the distribution function of maximum brightness over the decline rate at 
different z. 

What to do in this situation? First, the statistics must be improved to such an extent that, 
at a given high redshift, one could compare the distributions of SNe on the maximum luminosity 
- decline rate diagram. I.e., at each z one should build the distributions of the decline rates in a 
broad range of the peak luminosity. The hardest thing is, of course, to observe the intrinsically 
weak SNe la for sufficiently long time. The distributions found for high z could coincide with 
the one for nearby supernovae, but could also differ from it. The latter will mean that nearby 
statistical relations are not applicable to the distant events. But even in the case when the 
brightest supernovae at, for instance, z = 1 have the same distribution of their light curve slopes 
as the brightest samples at z < 0.1, this will tell no word about the coincidence of their absolute 
magnitudes if one does not use any additional data. 

Second, the individual differences are larger in the / band. This spectrum interval is more 
sensitive to the explosion physics. Perhaps, the influence of the burning mode on the / band 
will become more transparent in future theoretical investigations. The theorists should simulate 
the light curves in / for extended range of models and compare them with the observations. 
It is hard to measure the rest frame / flux for large z, since then it goes to the far infrared 
on Earth. But this must be possible with new orbiting telescopes, like the Next Generation 
Space Telescope (Dahlen & Fransson 1998, Rees 1998). The observations in / should help to 
put useful constraints on the supernova models for each real explosion. But the thorough study 
is still needed both to understand the role of the non-LTE effects in the emission of SNe la 
in near infrared and to develop the realistic 3D models of thermonuclear burning. Ultimately, 
the radiative transfer computations for those 3D-models will show the level of accuracy of ID 
models considered here. 
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Table 1: Parameters of the models 



model 


burning mode 


M WD (M ) 


-E51 


M 


36Ni (M ) 


opacity 


W7 


deflagration 


1.3775 


1.20 




0.60 


scattering 

+ absorption 


W7a 


deflagration 


1.3775 


1.20 




0.60 


pure absorption 


DD4 


delayed 
detonation 


1.3861 


1.23 




0.626 


scattering 

+ absorption 



Table 2: Chemical composition of the models (all masses are in units of M & ] 



56 Ni 

He 

C 

o 

Ne 

Mg 

Si 

S 

Ar 

Ca 

Fe 

58 Ni 



W7 (W7a) 
5.99 • 10- 1 
2.88 • 10~ 3 
4.65 • 10~ 2 

1.41 • 10- 1 

4.42 • 10~ 3 
8.64 • 10~ 3 
1.52 -10- 1 

8.62 • 10~ 2 
1.60 • 10~ 2 
1.24 • 10~ 2 
1.45 • 10" 1 

1.63 • IP' 1 



DD4 



6.26 • 10 _1 


3.28 • 10~ 3 

1.01 -lo- 1 

• 10~ 3 
■10~ 3 

■ 10- 1 

■ 10- 1 

•10~ 2 



1.20 

8.97 
2.58 
1.56 
3.47 



3.75 • 10~ 2 

1.59 -10- 1 





Table 3: Rise time to maximum light in different colors 



model 


rise 


M B 


rise 


M v 

-'"max 


f bol 

rise 


M bo1 

-'"max 


log^ax 


DD4 


17.1 


-18.91 


22.1 


-19.03 


9.4 


-19.36 


43.23 


W7 


16.0 


-18.71 


20.5 


-18.82 


8.6 


-19.38 


43.24 


W7a 


11.9 


-18.97 


13.9 


-19.18 


8.9 


-19.42 


43.25 
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Figure 1: The density and the velocity distributions on mass at the beginning of our calculations 
of models DD4 and W7. The density is plotted relative to the central value. We have started 
from p c = 9.76 • lCT 9 g cm~ 3 at t = 8.16 • 10 4 s for DD4, and from p c = 7.42 • 10~ 3 g cm~ 3 at 
t = 8 • 10 2 s for W7. 
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Figure 2: The distribution of the most abundant chemical elements throughout the ejecta for 
DD4 (top) and W7 (bottom). 
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Figure 3: The simulated UBVI light curves of DD4 (solid line) and W7 (dotted line). Asterisks 
and crosses correspond to the light curves of the SN 1998bu and SN 1992A, respectively. The 
observed light curves are shifted along the time axis so that their maxima in B occur at the 
same time as the maximum of the calculated flux from DD4. 
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Figure 4: The temperature distribution for DD4 (solid) and W7 (dashed) on 35th day after the 
explosion. The solid dots show the locations where the optical depth measured from the surface 
reaches the value of 2/3. 
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Figure 5: The comparison of the light curves for two runs of the model W7. a) The extinction 
is treated as pure absorption and thereby the expansion line opacity is fully (and incorrectly) 
included into the energy equation (the model W7a, dotted). b) Scattering is taken into 
account; the expansion effect in line opacity is not included in the energy equation, i.e. the 
correct equation is solved (the model W7, solid.) 
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